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Abstract 

Available analyses on moments of the spin structure function <?i use different methods and are barely 
consistent with each other. We present a higher twist analysis of using a method consistent with 
the studies of T™ and Tj~ n already published. The twist-4 coefficient /2 is extracted. One result is 
that the higher twist coefficients seem to alternate signs: the relatively larger twist-6 contribution is 
partly suppressed by the twist-4 and twist-8 contributions. The size of twist-6 can be due to the elastic 
contribution to the moments. 

High precision data on doubly polarized electron-nucleon scattering from Jefferson Lab (JLab) have been 
analyzed in the transition regime from asymptotically free to strongly interacting quarks |1I2I3I4| , Studying 
quark-gluon and quark-quark interactions is important to understanding quark confinement. Such study can 
be cast in the Operator Product Expansion (OPE) formalism, which describes in particular the evolution of 
structure functions and their moments. The Cornwall-Norton moment is the integral of the structure function 
over x. Here, x = Q 2 /2Mv is the Bjorken variable, Q 2 is the four-momentum transfer from the electron to 
the nucleon, v is the energy transfer and M is the nucleon mass. In OPE, the first moment of g\{x, Q 2 ) can 
be written as: 



r 1 (g 2 ) s / 1 dx 3l (x,Q 2 ) = E r=2 ,4.. 



MQ 2 ) 



where the n T (Q 2 ) are sums of twist elements added up to twist r. The twist is defined as the mass dimension 
minus the spin of an operator. Twist elements> 3 can be related to quark-quark and quark-gluon interactions. 
Hence they are important quantities for confinement study. The leading twist coefficient is: 

nl {n) {Q 2 ) = C ns (Q) 2 (ia 8 ± ±g A ) + C s {Q) 2 \a Q 

where C ns and C s are flavor non-singlet and singlet Wilson coefficients that represent the Q 2 -dependence due 
to QCD radiations [U, g a = 1.267(35) is the triplet axial charge a§ ~ 0.579(25) is the octet axial charge 
Hi] and a is the singlet axial charge. In the MS renormalization scheme that will be used here, ao = AE 
where AS is the contribution of the quarks to the nucleon spin. The next to leading order twist coefficient 
is: 



MQ 2 ) = MQ2) + Ad 2 (Q 2 ) + 4/ 2 (Q 2 )) 

a 2 (cfo) is a twist two (three) target mass correction that can be related to higher moments of g\ (of g\ and 
52), and fi is the twist four contribution [Jj. 

OPE analysis can also be carried out using Nachtmann moments [Hj, in which the target mass corrections 
are done by an appropriate combination of g\ and gi in the moment's definition. Such an analysis of 
the new JLab EGla data have been carried out on Tj ^J. On the other hand, analysis of T™ ^H] and 
the flavor non-singlet T\~ n pi] were done using Cornwall-Norton moments. The results at Q 2 = 1 GeV 2 
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for h are / 2 P = 0.039 ± 0.022(stat)i , ; 'i8(syst)±0.030(low x)Zo™[(a s ), fi = 0.034 ± 0.005 ± 0.043 and 
fP- n = -0.13 ± 0.15(uncor.)t l ; l ^(cor.) where uncor. (cor.) specifies the error due to the uncorrelated 
(correlated) experimental uncertainty. The /jq results are /ig * /Af 4 =0.011±0.013(stat)l°;°J°(syst)±0.011(low 
x)±0.000(a s ), (J%/M 4 = -0.019 ± 0.002 ± 0.017 and $~ n /M 4 = 0.09 ± 0.06(uncor.)±0.01(cor.) where the 
asterisk in ^g* recalls that this coefficient contains only a twist 6 term. 

These results, while coming from the same set of data, barely agree. The disagreement could come from 
the fact that the low-a; extrapolation procedures differ in the three analyses, or the lower Q 2 considered for the 
fits are different (1 GeV 2 for p, 0.5 for n and 0.8 for p-n), or the target mass corrections are treated differently 
in the Nachtmann and Cornwall-Norton analyses: in the Nachtmann moments, target mass corrections are 
added to all orders while in the Cornwall-Norton analyses, only the first order is corrected for. The Cornwall- 
Norton analyses indicate that twist 4 and twist 6 terms are of similar magnitude (although twist 6 is larger) 
but opposite sign, leading to a partial cancellation of higher twist effects. This is not as clear from the 
Nachtmann analysis. To clarify this issue, it would be beneficial to provide consistent OPE analysis of the 
data. In that light, we have redone a Cornwall-Norton analysis of the T± data consistent with the r™ and 
r^~™ analyses. 

The low— x extrapolation of the JLab and world data was redone, as in the F™ and Y\~ n analyses, using the 
Thomas-Bianchi parametrization up to the invariant mass squared W 2 = 1000 GeV 2 . The uncertainty 
was estimated by varying all the parameters within their range given in ^21- A R- e gg e form was used 
beyond W 2 = 1000 GeV 2 on which an uncertainty of 100% was assumed. The elastic contribution to the 
moments was estimated using the parametrization of Mergell et al. 22 • A 2% uncertainty was assumed. The 
JLab EGla experiment (that will mainly determine the higher twist magnitude) is dominated by systematic 
uncertainties. Its point to point uncorrelated systematic uncertainties were separated from its correlated 
ones, and added in quadrature to its statistical uncertainty. This error was used in the OPE fit. The effect 
of the point to point correlated uncertainty was accounted for by shifting the EGla data set and using it as 
a new input in the fit. 

Fitting the world data for Q 2 > 5 GeV 2 and assuming no higher twist effects above Q 2 — 5 GeV 2 yields 
AS = 0.154 ± 0.066. The target mass correction ao(Q 2 ) — L dx [x 2 g\{x 1 Q 2 )) , where gi(x,Q 2 ) contains 
only a twist-2 contribution, was estimated with the parton distribution parameterization of J. Bluemlein 
and H. Boettcher J^j. The twist-3 contribution d 2 {Q 2 ) was obtain from the SLAC E155x experiment 
|16j . Although accounting for the Q 2 dependence had little effect on the fit, a Q 2 — dependence of the 

form A(Q 2 ) = A{Q 2 ) ) (a s (Ql) / a s (Q 2 )) was assumed for a (Q 2 ) and d 2 (Q 2 ) with b = -0.2 and b = -1 
respectively. Aqcd — 0.371q'q7 was used in computing a s (Q 2 ). 

The world data together with the OPE leading twist evolution (LT) of T\(Q 2 ) and the elastic contribution 
to T\{Q 2 ) are shown in the figure below. The band at zero is the point to point correlated uncertainty on 
the JLab EGla data. The dot-dashed line is the result of fit 1 (see table). 
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Figure 1: World data on T\{Q 2 ). The gray band (LT) is the pQCD leading twist evolution. The band on the 
horizontal axis is the point to point correlated uncertainty for the JLab CLAS experiment. The uncorrelated 
uncertainty is of the size of the square symbols. The error bars on the open symbols are systematic and 
statistic added in quadrature. The dash-dotted line is a fit of the data starting at Q 2 nin = 0.1 GeV 2 . 

To check the convergence of the OPE series, the lowest Q 2 value, Q^ in was varied, as well as OPE series 
truncated to twist-8. The results are given in the table below. All the twist coefficient values are given 
for Q 2 = l GeV 2 . The first error represents the uncorrelated uncertainty, coming mainly from the statistical 
uncertainty, and the second is the point to point correlated uncertainty. 



fit 


Q 2 . 

^ mm 


h 


M4/M 2 


im/m 4 




1 


1.0 


-0.138±0.024t°;^j 


-O.OS5±0.011±°;°™ 


0.110±0.014±°;°g 




2 


0.8 


-0.120±0.017+°oig 


-0.047±0.073+°;°Q7 


0.099±0.008+°;°32 




3 


0.8 


-0.144±0.057to'.?27 


-0.057±0.025to.o28 


0.124±0.058to°3? 


-0.014±0.032±g'°26 


4 


0.6 


-0.160±0.027+°;^g 


-0.064±0.012+°;g^ 


0.143±0.021+°;°^ 


-0.026±0.008+g;°}g 



All the fit results are very consistent with each others. In fits 4 and 6, the smallness of fjL% tends to indicate 
the convergence of the OPE series. There is good agreement between the T™ and Y\~ n analyses and our 
analysis, although the central values differ noticeably. Also, our results show the same trend as the results 
from the neutron ^01 an d Bjorken sum analysis [llj : The f% coefficient tends to display an opposite sign as 
the [Iq coefficient. The alternation of signs seems to continue with /is, which indicates that the overall effects 
of higher twist are suppressed. This would indicate that we should expect hadron-parton duality JHj to hold 
for g\ at the scale at which the higher twist coefficients were extracted. The fact that stands out as the 
largest coefficient is, candidly, not surprising since in our Q 2 ranges, the Q 2 -behavior is dominated by the 
elastic contribution which roughly behaves as l/Q 4 . This feature was also seen in the Bjorken sum analysis 
but not in the neutron analysis in which the elastic contribution is smaller. 

These results can be compared to non-perturbative model predictions: fa — — 0.037±0.006 [17) . ^a/M 2 = 
-0.040 ± 0.023 (QCD sum rules [201), f 2 = -0.10 ± 0.05 (MIT bag model [HI) and f 2 = -0.046 (instanton 
model (22l)- As for the extracted ji and /14, all the predictions are negative. The MIT bag model and QCD 
sum rules agree best with the fit results, although the other predictions are not ruled out. 

Although agreeing well within uncertainties, Y\ ^ T\~ n + T™ from the analyses ^I] and This comes 
from the fact that the AS extracted from the proton and neutron analysis are very different: AE^™* 1 = 
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0.15(0.35). This implies that the asymptotic values for the , F™ and rf " that anchor the OPE evo- 
lutions used in the fit are inconsistent. As an example, an offset of the Bjorken sum asymptotic value of 
(AS? - AE™ = -0.2) /9 changes the value of /f at Q 2 = 1 GeV 2 by a factor 2 and the value of [i\~ n by 
50%. 

From the result of fit 1, we can extract the proton color polarizabilities which are the responses of the 
color magnetic and electric fields to the spin of the proton |17I18| : x V e = —0.08 ± 0.02 i !os an d Xb = 
0.06 ± 0.08±^;o|. As for the neutron, these are of opposite sign and compatible with zero. 

The fact that the higher twist effects are small (at Q 2 = 1 GeV 2 ) is somewhat surprising and exciting. 
However, it implies that accurate measurements are more delicate: the size of the uncertainty is presently of 
the size of the central value itself. In particular, the high energy missing part, very substantial for the JLab 
data, introduces a significant uncertainty. The 12 GeV upgrade of JLab will improve on this issue and push 
the Q 2 -coverage of the measurements. It will help in measuring the higher twist coefficients precisely. 
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